Novel Techniques for Decomposing Diffuse Backgrounds 



Brandon S. Hensley/'B Vasiliki Pavlidou,^'[l| and Jennifer M. Siegal-Gaskins^-'^'llI 

'Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA 
^Max-Planck Institute for Radioastronomy, 53121 Bonn, Germany 
^California Institute of Technology, Pasadena, CA 91125, USA 
'* Center for Cosmology and Astro- Particle Physics, 
Ohio State University, Columbus, OH 43210, USA 
(Dated: October 30, 2012) 

The total anisotropy of a diffuse background composed of two or more sources, such as the Fermi- 
LAT-measured gamma-ray background, is set by the anisotropy of each source population and the 
contribution of each population to the total intensity. The total anisotropy as a function of energy 
(the anisotropy energy spectrum) will modulate as the relative contributions of the sources change, 
implying that the anisotropy energy spectrum also encodes the intensity spectrum of each source 
class. We develop techniques, applicable to any such diffuse background, for unraveling the intensity 
spectrum of each component source population given a measurement of the total intensity spectrum 
and the total anisotropy energy spectrum, without introducing a priori assumptions about the 
spectra of the source classes. We demonstrate the potential of these methods by applying them to 
example scenarios for the composition of the Fermi-LAT gamma-ray background consistent with 
current data and feasible within 10 years of observation. 

PACS numbers: 95.35.-l-d, 98.70.Vc, 98.54.Cm, 98.70.Rz 



I. INTRODUCTION 

Diffuse emission, from radio to gamma-ray frequen- 
cies, encodes a wealth of information about fundamen- 
tal physics, cosmology, and a variety of astrophysical 
systems. Prominent examples include the cosmic mi- 
crowave background, a snapshot of the very early uni- 
verse at microwave frequencies [l| ; direct and reprocessed 
starlight between infrared and ultraviolet wavelengths, a 
record of the star formation history of the universe |2s- 
3]; thermal emission from accretion processes in X-rays, 
which traces the growth of black holes through cosmic 
time Q ; and non-thermal emission in gamma rays from, 
e.g., blazars star- forming galaxies [9], millisecond 

pulsars llOl, an d, possibly, annihilating or decaying dark 
matter [llFlil. 

Both the intensity spectrum and the degree of 
anisotropy of the diffuse signal have been successfully 
used to uncover valuable information about the physics 
and astrophysics of the processes and sources that are 
responsible for the diffuse emission in each waveband, 
the cosmic microwave background bein g th e most cel- 
ebrated example of both methods [e.g.. IT^. [TBj . How- 
ever, when more than one source class contributes to the 
diffuse emission, complications to such analyses arise. 
Traditionally, determining the individual contributions 
of source classes in a multi-population diffuse signal has 
relied on careful modeling and subtraction of intensity 
spectra [l6j . This process is hindered by systematic un- 
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certainties in the theoretical understanding of astrophysi- 
cal source classes, which limit our ability to detect a pos- 
sibly subdominant signal, especially one with unknown 
or poorly constrained properties. 

In this work, we illustrate that a measurement of the 
anisotropy of a diffuse background at a fixed angular 
scale as a function of energy, when combined with the 
total intensity spectrum of the background, can be used 
to decouple the contributions to the background of each 
source population, thereby yielding an intensity spec- 
trum for each source class. We focus here on applica- 
tions to the isotropic diffuse gamma-ray background, as 
its large range of plausible compositions allows us to il- 
lustrate many of the techniques presented here. 

The isotropic gamma-ray background (IGRB), the dif- 
fuse gamma-ray emission at energies above ~ 100 MeV 
that is isotropic on large angular scales, is one of the most 
promising observational targets for the discovery of new 
physics in this decade, such as a signature from dark 
matter annihilation or decay. Alt houg h the IGRB has 
been observed since the 1970s [13, llSj the Large Area 
Telescope (LAT) aboard the currently operational Fermi 
Gamma-Ray Space telescope {Fermi 'l9'|) is improving 
both the energy range and the angular accuracy of these 
observations. The LAT also resolves more bright point 
sources than previous missions due to its increased sensi- 
tivity, providing valuable information about gamma-ray 
source populations via detected members. As a result, 
the LAT collaboration has reported a more precise mea- 
surement of the IGRB intensity spectrum [20[, and for 
the first tinie has measured the small-scale anisotropy of 
the IGRB glj. 

In addition to any possible exotic signal, confirmed as- 
trophysical gamma-ray sources such as gamma-ray loud 
active galactic nuclei (blazars) and star-forming galaxies 
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are guaranteed to contribute significantly to the IGRB at 
some energy. Reference [2^ showed that by combining 
the spectral and anisotropy properties of the IGRB, it is 
possible to identify the presence of a second, even sub- 
dominant, component, such as a signal from dark mat- 
ter annihilation or decay, over a dominant, astrophysi- 
cal contribution. We extend this approach by developing 
techniques that allow the intensity spectra of the individ- 
ual components to be reconstructed without requiring a 
model or prediction for any of the contributions. 

For diffuse backgrounds composed of emission from un- 
correlated source populations, we show that under cer- 
tain conditions, if the intensity energy spectrum (differ- 
ential photon intensity as a function of energy) and the 
anisotropy energy spectrum (angular power at a fixed 
multipole as a function of energy) of the diffuse back- 
ground are both measured with sufficient accuracy, the 
shape of the intensity energy spectrum of each compo- 
nent can be recovered; in some cases the absolute nor- 
malizations of the intensity spectra are also recoverable. 
Similarly, in some cases the amplitude of the angular 
power spectra of the individual components can also be 
determined. We discuss the conditions under which such 
decompositions are feasible, and demonstrate these novel 
techniques on plausible scenarios for the IGRB composi- 
tion. Although our examples are restricted to the IGRB, 
the methods presented here can be applied to any diffuse 
background at any wavelength. 

In S|ll] we introduce the formalism common to all of 
the decomposition techniques. In mil] we define IGRB 
component models and simulated observations used for 
the example scenarios. The details of each decomposition 
technique are described in i jlVI example scenarios illus- 
trating a subset of the techniques are also presented. We 
extend our approach to selected three-component scenar- 
ios in fjVl We discuss the potential of these techniques for 
understanding gamma-ray source populations in WII 



II. TWO-COMPONENT DECOMPOSITION: 
METHODS 

The two properties of diffuse emission we will use are 
the differential intensity energy spectrum I{E) (photons 
per area per time per solid angle per energy) and the 
angular power spectrum Ci of a sky map of the inten- 
sity. The angular power spectrum is defined as — 
(|a^mP), where aim are the coefficients of the expansion 
of the intensity map in the basis of spherical harmonics. 
We also define the fluctuation angular power spectrum 
Ce = Ce/P, where / is the mean intensity of the emis- 
sion with intensity angular power spectrum Cg. Because 
Ci describes fluctuations in units of the mean, the fluc- 
tuation angular power at a fixed £ is energy-independent 
for a signal arising from a single population of sources 
with identical observer-frame intensity spectra. In the 
following we assume that each distinct component of the 
diffuse emission meets this criterion. 



Variation between the source spectra of individual 
members of a population can result in fluctuation an- 
gular power which is energy dependent because the rel- 
ative contributions of spectrally different sources within 
a population change with energy (e.g., harder sources 
contribute relatively more flux at high energies than at 
low energies). In addition, for cosmological source popu- 
lations, energy-dependent fluctuation angular power can 
also arise due to redshifting of sharp features in the source 
spectra, such as line emission or abrupt cut-offs (see, 

e.g., mM)- 

In practice, if a component of the emission arises from 
a population of sources, we assume that the requirement 
that the single-population Ci is energy-independent is 
satisfled if the variation in the intensity spectra of in- 
dividual members of the population is sufficiently small 
that the deviation of the fluctuation angular power from 
an energy-independent quantity is at a level smaller than 
the uncertainty on the anisotropy measured by a speciflc 
observation. We comment on the validity of this assump- 
tion in the context of the IGRB in mil[ 

Our approach exploits the energy independence of the 
single-component fluctuation angular power, and so it is 
convenient for us to work with Ci. If we consider a sce- 
nario in which the diffuse emission is composed of emis- 
sion from two spatially uncorrelated components with in- 
tensity spectra Ii (E) and I2 (E) and angular power spec- 
tra Ci^i and Ci^2, then the total intensity is simply the 
sum of the two components, 

Itot{E) = hiE) + hiE) . (1) 

The angular power spectrum of the total signal for un- 
correlated components is the sum of the angular power 
spectra of the components, 

Ci,tot{E)^Ce^i{E)+Ce^2{E). (2) 

Rewritten in terms of the fluctuation angular power. 



Ci^tot{E) 



h[E) 
hot{E) 



h{E) 
Itoi{E) 
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(3) 



This is the fluctuation anisotropy energy spectrum for 
the case we consider. In the following we will always 
use the term "anisotropy energy spectrum" to refer to 
the fluctuation angular power of the total emission as a 
function of energy. 

With sufRcient photon statistics, /tot and Cg^tot can be 
determined at each energy from observations. If there is 
a way to also determine Cg^i and Ci^2 from the data, we 
can solve Eqs. ([T|) and for Ii and I2: 



1 1 — hot 



Ci^2 ± \ Ci^iCi^tot + CiaCixot — Ci^iCi', 



Ci.i -f Ci 



1,2 



Ci^i T v Ci^iCixot + Ci^2Ci^tot — Ci^iCi' 

CiA + Cl2 



(4) 



(5) 
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TABLE I: Summary of two-component decomposition techniques. 



Method 


Observational Signature 


Inferred Properties of 
Components 


Intensity 

Normalization 

Recovered? 


Fluctuation 
Angular Power 
Recovered? 


Double plateau 


Plateaus at both high and 
low energies observed in 
anisotropy energy spectrum 


One source dominant in 
anisotropy at low energies, 
other source dominant at 
high energies 


Yes 


Yes 


Low-Anisotropy 
Plateau 


Anisotropy energy spectrum 
rises from (falls to) a low- 
anisotropy plateau at low 
(high) energy 


Source that is subdominant 
in intensity is much more 
anisotropic than the domi- 
nant source 


No 


No 


High-Anisotropy 
Plateau 


Anisotropy energy spectrum 
falls from (rises to) a high- 
anisotropy plateau at low 
(high) energy 


Source that is subdominant 
in intensity is much less 
anisotropic than the domi- 
nant source 


Yes 


No 


Known 

Zero-Anisotropy 
Component 


None; requires a pnori 
knowledge that one of the 
two components is isotropic 


One source is completely 
isotropic 


No 


No 


Minimum 


Minimum observed in the 
anisotropy energy spectrum 


Both source components 
have comparable intensity 
and anisotropy such that 
Eq. [2S1 is satisfied at some 
energy 


Yes 


Yes 


Multiple-^ 
Measurements 


Two distinct anisotropy en- 
ergy spectra can be obtained 
at two different £ 


Ce is a function of £ for at 
least one source such that 
two distinct anisotropy en- 


Yes 


Yes 



ergy spectra can be obtained 
at different £ 



If there is an energy ~ Eq around which only one com- 
ponent is expected to contribute to the total intensity 
(i.e., an energy range around Eq where hiEo) / Itot{Eo) « 
0), the anisotropy energy spectrum will be flat over this 
energy range. Then from Eq. ([3]) we immediately obtain 
Ce,i — Ce,tot{Eo) from the anisotropy of this baseline. 
A similar flat baseline could result if two source classes 
have the same spectral shape over an energy range, but 
such a scenario is unlikely for the source classes consid- 
ered here. In each of the following cases, we will assume 
either that we can obtain the Ci of one of the two source 
classes in this way, or that one source class is known to 
have Ce — 0. 

We discuss six distinct two-component decomposition 
techniques below. Some of them allow us to extract the 
component intensity spectra, while others only allow us 
to derive the shapes of one or both intensity spectra up 
to unknown normalization constants. Some of the tech- 
niques also yield measurements of the fluctuation angu- 
lar power spectra of each component source population. 
Table U gives a summary of these techniques and their 
applicability conditions. We emphasize that in all cases 
we make the following three assumptions: (1) the diffuse 
background is composed of emission from uncorrelated 



source classes, (2) the fluctuation angular power of each 
individual component is independent of energy, and (3) 
the fluctuation angular power of one component can be 
directly measured from the data at some energy or is 
known to be negligibly small. 

III. PARAMETERS OF EXAMPLE SCENARIOS 

A. IGRB Component Models 

To illustrate the decomposition techniques, we apply 
them to example scenarios that could be measured by 
the Fermi-LAT within 10 years of observation time. Each 
scenario is a two-component IGRB model consistent with 
current observations of the measured IGRB intensity en- 
ergy spectrum 2Q] and the measurement of the IGRB 
anisotropy energy spectrum [2l|. In particular, we con- 
sider a power law component with slope and fluctuation 
anisotropy that describes the IGRB well at low energies 
but that may break at high energies, as well as a Galac- 
tic dark matter annihilation component with one of two 
benchmark spectra (annihilation to a t+t^ or bb final 
state). 
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The power law component encapsulates likely contri- 
butions from several source classes such as blazars, but 
we assume, in accord with the data, that this emission 
can be characterized by a single Ci. The adopted Ci 
in this work differs from that derived in [2l[ as we use 
the intensities reported by [20| rather than those used 
in the anisotropy analysis, which were subjected to less 
cleaning. As this cleaning likely removed only isotropic 
contamination, we expect that the angular power Cp de- 
rived in [U also describes the IGRB of [13] , and thus we 
check for consistency with Cp only. 

Although we do not consider them explicitly in our ex- 
ample scenarios, other known gamma-ray source popula- 
tions which may contribute significantly to the intensity 
and/or anisotropy of the IGRB at some energies include 
star-fori ning galaxies d, [25l - [29| . gamma-ray loud radio 
galaxies |3C| and Galactic millisecond pulsars [13, HJ . 

For source classes relevant for the IGRB, the three as- 
sumptions stated above, which are necessary to imple- 
ment our methods, are very likely to be valid. While all 
cosmological populations will generally trace large-scale 
structure, on the small angular scales considered in the 
Fermi LAT anisotropy analysis [2l| we do not expect 
strong spatial correlations between, e.g., blazars and ex- 
tragalactic dark matter ^32], and naturally the angular 
distribution of any Galactic source population is entirely 
uncorrelated with that of any extragalactic population. 
For these source classes, the variation between observer- 
frame source spectra within the population is likely to 
be sufficiently small (or zero, in the case of Galactic dark 
matter annihilation or decay) to induce at most a mild 
energy dependence in the fiuctuation angular power, al- 
though we caution that a careful investigation of the ex- 
pected magnitude of this effect is needed. Furthermore, 
most expected cosmological contributors to the IGRB are 
not expected to exhibit features in their spectra which 
are sharp enough to lead to significant energy depen- 
dence of the fluctuation angular power due to redshift- 
ing. Although in some dark matter models sharp fea- 
tures are present in the photon spectra, in many scenar- 
ios the Galactic dark matter signal is expected to dom- 
inate over the extragalactic signal in both intensity and 
anisotropy, so any energy dependence in the fluctuation 
angular power of the extragalactic dark matter compo- 
nent would likely have a subdominant effect. Finally, the 
expected contributors to the IGRB have different spec- 
tral shapes, and consequently it is not unlikely that at 
certain energies all but one component will supply a neg- 
ligible contribution, in which case the fluctuation angular 
power of the dominant component could be measured as 
described in the previous section. 

The level of the blazar contribution to the IGRB is un- 
certain, with different calculations spanning a large range 
of possibilities (e.g., (6l-l8l, l33l - l35j ). The strongest bounds 
on the contribution of blazars to the intensity of the 
IGRB have been obtained by requiring that the adopted 
model for the blazar population does not exceed the mea- 
sured IGRB anisotropy [§3, HH] ; these bounds limit the 



blazar IGRB intensity contribution to < 20% in the 1- 
10 GeV band. The spectral shape of the blazar contribu- 
tion is dependent primarily on the distribution of blazar 
spectral indices in the gamma-ray range [SG^. Since BL 
Lac-type blazars and flat-spectrum radio quasars (FS- 
RQs), the two largest subclasses of blazars, generally 
have different spectral properties, the shape of the col- 
lective intensity spectrum depends on the relative abun- 
dances of these two subclasses in the unresolved blazar 
population, which is uncertain due to the difficulties in 
obtaining BL Lac redshifts and assessing the prevalence 
of BL Lacs in the high-redshift universe [3] . Additional 
uncertainties enter through considerations regarding the 
fraction of blazars with spectral breaks [s^] ■ 

Source intensity spectra at energies above a few tens 
of GeV from high-redshift populations are attenuated by 
interactions with the EBL, which consists of infrared, op- 
tical, and ultraviolet photons primarily from direct and 
reprocessed starlight throughout cosmic history. While 
Fermi observations have produced constraints on models 
of the EBL [s^ [s^], the details remain quite uncertain 
[1, 140I - I42I . In some of the models we consider, we use 
a broken power-law model for the non-dark-matter emis- 
sion, which can act as a proxy for EBL attenuation as 
well as for accounting for the intrinsic properties of the 
emitting sources. 

We assume an anisotropy for our composite power law 
consistent with observations rather than tying it to mod- 
els of a given source class. A significant contribution 
to the anisotropy is expected to come from unresolved 
blazars. The anisotropy properties of the blazar contri- 
bution to the background are generally dependent on the 
details of the blazar luminosity function, and therefore 
similarly uncertain as the overall amplitude of the col- 
lective blazar intensity. The derived anisotropy of the 
IGRB is comparable to the level of blazar fiuctuation 
anisotropy predicted by theoretical work (see, e.g., ,32^]), 
though blazars are expected to be too anisotropic to con- 
tribute the entire IGRB [Mllil. 

For the dark matter component in our example sce- 
narios, we model the emission from pair annihilation of 
WIMP dark matter particles in Galactic subhalos. We 
consider the photon intensity spectra given in [i^l pro- 
duced by annihilation into two benchmark final states: 
(1) 66, which generates a relatively soft continuum photon 
spectrum primarily from the decay of neutral pions pro- 
duced by the hadronization of quark jets, and (2) r+r~, 
which produces a harder photon spectrum due to a sig- 
nificant contribution from final state radiation associated 
with the production of charged leptons. We choose values 
of the annihilation cross section between the canonical 
value for a thermal relic {cfv)q = 3.0 x lO^^cm'^s"^ [il] 
(see also [45j) and 33 times that value. Different con- 
straints on dark matter models can be obtained under 
different assumptions for various targets, (e.g., (icl - fs^ ). 
The fluctuation anisotropy from dark matter annihila- 
tion, which is determined exclusively by the spatial distri- 
bution of the dark matter, has been predicted for annihi- 
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lation in Galactic dark matter subhalos by Refs. [53l - l56f . 
We present three models with a dark matter component: 
one that falls within the typical predictions for the dark 
matter anisotropy and two models that do not. How- 
ever, we emphasize that the model parameters adopted 
for each example scenario were chosen to be illustrative of 
the decomposition methods rather than to represent the 
most plausible compositions of the i^ermi-LAT IGRB. 



B. Error Analysis 

We compute error bars for the example IGRB inten- 
sity and anisotropy energy spectra assuming observations 
with the Fermi-hAT. The la error bars for the total 
fluctuation anisotropy in each energy bin were computed 
using the formula [57|: 



{2i + 1) AiUy 



(6) 



where d is the total fluctuation angular power spectrum, 
A£ is the width of the multipole bin, /sky is the fraction 
of the sky used to calculate the angular power spectrum, 
C'n = (47r/sky/-/V.^) is the power spectrum of the photon 
noise associated with the total measured emission, with 
iV^ the total number of photons collected during the ob- 
servation period in the sky region analyzed, and Wi is 
the beam window function of the instrument. 

We approximate the PSF of the LAT as a circular 
Gaussian beam with energy-dependent width (7b (E) , de- 
termined from the 68% containment angle radius re- 



ported in the P7_V6 performance curves [61|, so We 



exp (—£^(7^/2), which is the window function of a Gaus- 
sian beam of width ah- We evaluate ah at the log cen- 
ter of the energy bin. For the example scenarios, we 
show anisotropy energy spectra at £ — 175, and take 
A£ — 50, choices made to ease comparison with the re- 
sults reported in the Fermi anisotropy analysis f^T] which 
used A£ = 50 and focused on the Poisson angular power 
measured at £ > 150 to limit contamination from Galac- 
tic diffuse emission. 

Following the Fermi anisotropy analysis, we assume a 
sky fraction /sky = 0.32 is used to perform the anisotropy 
and intensity measurements, i.e., a large fraction of the 
sky is masked. We take the field of view of the LAT to 
be = 2.4 sr, and approximate the energy-dependent 
effective area of the LAT from the reported performance 
curves. For a specified all-sky observation time iobsj we 
calculate the number of photons detected outside the 
mask to be Nj = J^^"^ dE ^ Acg{E) /sky tobs, where 

■j^ is the total (energy-dependent) differential intensity 
of the IGRB, A^s is the energy-dependent effective area 
of the instrument, and we have assumed observations in 
all-sky survey mode and uniform sky exposure. 

The errors on the intensity energy spectra represent 
the Poisson noise associated with the number of photons 



collected in each energy bin, as well as an assumed 20% 
uncertainty on the effective area of the instrument as es- 
timated from the performance curves; we note that the 
uncertainty on the effective area does not contribute to 
the uncertainty in the fluctuation angular power spectra 
since normalizing intensity fluctuations to the mean map 
intensity removes the effective area from the calculation. 
For the all-sky observation time of 10 years assumed in 
the examples presented in this work, we find that for the 
adopted energy binning, the uncertainty in Agg domi- 
nates the error bars on the intensity spectra up to a few 
hundred GeV for all of the scenarios considered here, re- 
sulting in relatively uniform error bars as a function of 
energy. 

The errors on the simulated intensity and anisotropy 
spectra can be propagated through the decomposition 
equations, enabling calculation of error bars for the 
decomposed intensity spectra. Often, the Ce of one 
source class must be obtained from a fiat baseline in the 
anisotropy energy spectrum over the energies where that 
source class dominates the anisotropy. An estimate for 
the Ce of such a source class can be obtained by taking 
the weighted mean of the baseline points with error equal 
to the weighted error of the mean. The decomposed in- 
tensity spectra can then be written in terms of quantities 
for which a mean and error bar can be computed. The 
la and 3a confidence intervals are determined via Monte 
Garlo. 



IV. TWO-COMPONENT DECOMPOSITION: 
APPLICATIONS 

The formalism for each of the six two-component de- 
composition techniques summarized in Table U is given in 
the following subsections. 



A. Double Plateau 

As previously described, we can infer the value of Ce.i 
by observing a flat baseline (plateau) in the anisotropy 
energy spectrum at either low or high energies. In the 
event that we observe a second such plateau in the 
anisotropy energy spectrum (Fig. [T]) with amplitude ei- 
ther above or below the level of the first plateau Ce.i {i.e. 
where Ii <C I2), we can obtain a value for C'e,2. In this 
case Eqs.|l]and[n]for the intensity spectra of the two com- 
ponents can be solved directly. This is a double plateau 
decomposition, and corresponds to the case that one of 
the components dominates the anisotropy at low ener- 
gies, while the other dominates the anisotropy at high 
energies. In this case we assume the anisotropy energy 
spectrum increases or decreases monotonically between 
the two plateaus; the case of local extrema is discussed 
in the case of a minimum decomposition in ijlVEl A dou- 
ble plateau is a particularly ideal case because the input 
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spectra can be derived exactly, without making any as- 
sumptions about the relative contributions of the source 
classes to either the total intensity or anisotropy that 
cannot be inferred directly from the observed spectra. 

As an example scenario, shown in Fig. [1] we choose 
™dm = 300 GeV, (crw) = 20{av)o, and annihilation into 
T+r^ for the dark matter intensity spectrum, and adopt 
a broken power law for the remaining intensity spectrum. 
The anisotropics were taken to he Ce — 1.2 x 10^^ sr, 
and Ci = 5 X 10~^ sr for the broken power law and dark 
matter signals, respectively. Because the anisotropy en- 
ergy spectrum is still rising slightly between the last two 
data points, the estimate for the of the dark matter 
component will be biased low, thus slightly biasing the 
decomposed spectra away from the true value. 



B. Low- Anisotropy Plateau 

We now consider a scenario in which a low-anisotropy 
plateau is measured in the anisotropy energy spectrum 
at low or high energies, and the anisotropy rises from the 
plateau at low energies or falls to the plateau at high 
energies. This corresponds to a case where the compo- 
nent subdominant in intensity at the plateau has a much 
higher anisotropy. 



with A > 1. 

In this case, Eq. ^ can be written as 



(7) 



e,tot 



= 1- 



^tot 



h_ 

hot 



or 



i,tot 



2^ + (l+A) (A 
hot \Uot 



AC,,! (8) 



(9) 



Now we define a; (i?) = h/hot anda;(£;) = CLtot/Ci,i~l, 
which can be determined by observations at each energy. 
Then we have 



with solution 



1± Vl + (1 + A)c 
1-l-A 



(10) 



(11) 



Since A 3> 1, as long as w > 1 we can approximate this 

by 



l±^(l + A)a 
1 + A 



Vl + A 



(12) 



where we have selected the + solution since a; is a non- 
negative quantity. Since lo is an observable, it is always 
possible to determine whether the a; > 1 condition holds. 



The shape of the subdominant spectrum can thus be de- 
rived up to a multiplicative constant. 

As an example scenario, shown in Fig. [21 we choose a 
dark matter particle with ttidm — 1000 GeV that annihi- 
lates to bb with {(jv) ~ 33(crw)o and a broken power law 
component to the intensity spectrum. We set the fluctu- 
ation angular power to = 1 x 10^^ sr for the broken 
power law component, and Ci = 2.5 x 10^'^ sr for dark 
matter. 



C. High- Anisotropy Plateau 

If a high-anisotropy plateau is measured in the 
anisotropy energy spectrum at low or high energies, and 
the anisotropy falls from the plateau at low energies 
or rises to the plateau at high energies, then a less 
anisotropic source must be making an increasing con- 
tribution to the background at energies far from the 
plateau. We now consider this scenario, corresponding 
to the case that one component is everywhere dominant 
in the intensity. 



and also more anisotropic 

In this case, Eq. [3]can be approximated by 
\hot J 

Immediately then we have 



t,tot 



and 



'a 



l,tot 



(13) 
(14) 

(15) 
(16) 
(17) 



The assumption that the higher anisotropy source is 
dominant must be satisfied at the energies where the 
high-anisotropy plateau is measured. Each subsequent 
point moving away from the plateau in energy will yield 
a value for the fractional contribution of the subdominant 
source via the ratio of Eqns. [16] and [171 Therefore, the 
appropriateness of the assumption that the anisotropic 
source is dominant can always be verified. This decom- 
position yields both the shape of the intensity spectra as 
well as their normalizations. 

As an example scenario, shown in Fig. [31 we choose a 
dark matter particle with ttidm — 200 GeV that anni- 
hilates to T+r~ with {av) = 6.7(fTu)o and a component 
with a power law intensity spectrum. We set the fluctu- 
ation angular power to (7^ = 1.2 x 10~^ sr for the power 
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FIG. 1: Example double plateau decomposition. Left: Total IGRB intensity (top), anisotropy energy spectrum (middle), and 
angular power Cp (bottom). Energy bins are equally spaced in log space, and all quantities are reported at the log center of 
the energy bin. Error bars assume tohs ~ 10 years of Fermi-LAT observations in sky-survey mode. If a data point is within 
3(T of zero, we place a 3cr upper limit bar in addition to the la error bars. Red triangles indicate the Fermi IGRB intensity 
and angular power measurements. The dark matter intensity spectrum corresponds to a moM ~ 300 GeV particle annihilating 
to r^r~. Right: The decomposed intensity energy spectrum of the power law component (blue circles) and dark matter 
annihilation (purple triangles) recovered using the plateau technique. The baseline points (open plot symbols) from which Ce,i 
and Ct^2 were determined were not decomposed. Each component's input intensity spectrum is overlaid in black. 



law component, and C'e = 5 x 10^^ sr for dark matter. 
These parameters are consistent with current observa- 
tional and theoretical constraints on the IGRB intensity, 
anisotropy, and dark matter properties. 

D. Known Zero- Anisotropy Component 

If one of the components of a two-component back- 
gromid is completely isotropic, then Eq. [3] becomes: 

C'^.tot = ( T~~ ) ^^.1 ^^^^ 
V-'tot / 

which may be rewitten simply as 

h = /totV^^W^- (19) 

Thus, in this case the shape of the spectrum of the com- 
ponent with nonzero anisotropy can be determined up 
to a multiplicative constant by measuring /tot(£^) and 
Ci^tot{E). This technique is a special case of the high- 
anisotropy plateau where Eqns. and [T7] hold exactly, 
regardless of which source is dominant. We emphasize 
that this method requires a priori knowledge that a zero- 
anisotropy component exists, but does not require knowl- 
edge of the shape or normalization of its intensity spec- 
trum. 



To demonstrate such a decomposition, as shown in 
Fig. |4l we analyze the intensity and anisotropy data 
as measured by Fermi-LAT; note that the measured 
anisotropy energy spectrum is consistent with no energy 
dependence. 

The decomposition presented in Fig. 2] yields spec- 
tra with slope consistent with the slope of the compos- 
ite IGRB. This decomposition technique recovers the 
shape but not the normalization of the intensity spec- 
trum of the anisotropic component, hence the normal- 
izations shown in Fig. 2] are arbitrary. We caution the 
reader that this is not a definitive decomposition of the 
Fermi-LAT IGRB due to uncertainty in matching the in- 
tensity and anisotropy measurements from two different 
studies, and in addition it is subject to the assumption 
that a zero-anisotropy component exists in the measured 
IGRB. 



E. Minimum 

Assuming Ci^i can be inferred from the data by mea- 
suring a plateau in the anisotropy energy spectrum, an 
additional way to determine Ce^2 and decompose the ob- 
served intensity spectrum exists if a local minimum is 
observed in the anisotropy energy spectrum. The condi- 
tion for observing a minimum in the anisotropy energy 
spectrum can be obtained by differentiating Eq. ^ with 
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Energy [GeV] Energy [GeV] 

FIG. 2: Example low-anisotropy plateau decomposition. Left: Total IGRB intensity (top), anisotropy energy spectrum (middle), 
and angular power Cp (bottom). Energy bins are equally spaced in log space, and all quantities are reported at the log center 
of the energy bin. Error bars assume tobs = 10 years of Fermi-LAT observations in sky-survey mode. If a data point is within 
3cr of zero, we place a 3a upper limit bar in addition to the la error bars. In this scenario the IGRB is composed of emission 
from a broken power law component and Galactic dark matter annihilation. The dark matter intensity spectrum corresponds 
to a TTiDM ~ 1000 GeV particle annihilating to bb. Right: The decomposed intensity energy spectra of the power law component 
(blue circles) and dark matter annihilation (purple triangles) are recovered using the low-anisotropy plateau technique. The 
baseline points at low energy (open plot symbols) from which Ci^i was determined were not decomposed. Each component's 
input intensity spectrum is overlaid in black. Red triangles indicate the Fermi IGRB intensity and angular power measurements. 
Note that this method recovers the normalizations of each component's intensity spectrum up to a multiplicative constant; the 
constant has been set to the true value to facilitate comparison with the input spectral shapes. 



respect to energy. We find that a minimum occurs when 

At the energy at which the local minimum occurs i?min, 
we can simultaneously solve this equation with Eq. [T] at 
the same energy to obtain 

^1,2 — ~ ~ • V^lj 

Substituting this relation into Eqs. S] and [5] completely 
determines the intensity energy spectrum for both source 
classes. In this case, a minimum decomposition is possi- 
ble. For appropriate levels of anisotropy and fractional 
contribution to the background, Eq. |3] states that the 
total anisotropy can be less than the anisotropy of ei- 
ther source class. Indeed, the minimum will exist only 
if there is an energy at which IiCi^i = l2Ci,2, and this 
requirement will always produce a minimum rather than 
a maximum in the total anisotropy since 




is non-negative. It also implies that a component which 
is always subdominant cannot be decomposed using 
this technique if its angular power is also smaller than 
that of the dominant component. A second caveat of 
this method is that a minimum can also occur when 
■j^ (irr) ~ which does not yield a constraint on Ci^2- 
This occurs when a low anisotropy component transitions 
from contributing an increasing fraction of the total in- 
tensity to contributing a decreasing fraction of the total 
intensity. Thus, the anisotropy of the resulting mini- 
mum must be intermediate between the anisotropy of 
each source class. On the other hand, the decomposition 
minimum occurs at an anisotropy below that of either 
source class. Hence, in practice, the two types of minima 
may be distinguished if the anisotropy energy spectrum 
is observed to take on values both above and below an 
observed baseline, in which case it must be a decomposi- 
tion minimum. 

The biggest challenge in the applicability of this tech- 
nique is that the uncertainties associated with an ob- 
served anisotropy energy spectrum will often be too high 
to allow a minimum to be measured. The depth of the 
minimum is given by the ratio of C^^tot at the minimum 
and the lesser of Ct^i and 0^,2- By rearranging Eq. [5T] 
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Energy [GeV] Energy [GeV] 

FIG. 3: Example high-anisotropy plateau decomposition. Left: Total IGRB intensity (top), anisotropy energy spectrum 
(middle), and angular power Cp (bottom). Energy bins are equally spaced in log space, and all quantities are reported at 
the log center of the energy bin. Error bars assume tobs = 10 years of Fermi-LAT observations in sky-survey mode. If a 
data point is within 3(7 of zero, we place a 3a upper limit bar in addition to the la error bars. In this scenario the IGRB is 
composed of emission from a power-law component and Galactic dark matter annihilation. The dark matter intensity spectrum 
corresponds to a mDM = 200 GeV particle annihilating to r"'"r~. Right: The decomposed intensity energy spectra of the 
power law component (blue circles) and dark matter annihilation (purple triangles) recovered using the low-anisotropy plateau 
technique. The baseline points at low energy (open plot symbols) from which Ce^\ was determined were not decomposed. Each 
component's input intensity spectrum is overlaid in black. Red triangles indicate the Fermi IGRB intensity and angular power 
measurements. 



and assuming Cg^i > Cg^2 we obtain 



^^4^ = „ (23) 

Thus the depth of the minimum varies between 0.5 
and 1 times Cf.i, with the largest depth occurring for 
Ci.i — Ct^2- For instance, although the example given 
in §IVBI (low-anisotropy plateau) in principle exhibits a 
local minimum in the anisotropy energy spectrum, the 
depth of this minimum would be very small and thus 
challenging to measure observationally. 

However, this does not mean that likely two- 
component scenarios for the IGRB that would allow a 
minimum decomposition do not exist. One such scenario 
can be envisioned if blazars typically exhibit spectral 
breaks. In combination with the different population- 
average spectral indices of BL Lacs and FSRQs, this 
could lead to a scenario in which over a certain energy 
range the IGRB is composed of emission from two source 
classes (BL Lacs and FSRQs), and the dominant con- 
tributor to the intensity transitions between the source 
classes at a few GeV (see, e.g., Fig. 3 of [131 )• If BL 
Lacs and FSRQs have approximately equal levels of fluc- 
tuation angular power (not unlikely, as Fermi has re- 
solved comparable numbers of sources in each of the two 



classes) , the minimum in the anisotropy energy spectrum 
would appear around the transition energy of a few GeV, 
where photon statistics would be fairly large and thus 
may allow a measurement of the minimum. 

We note that by examining the conditions under which 
a minimum occurs, we can also better understand the de- 
generacy in choice of sign in Eqs.|4]and[5] By rearranging 
Eq. El which gives 

hCi,2~hCt^^ (24) 

we see that the proper sign is determined based upon 
the sign of l2Ci,2 — IiCg^i- As the relative contributions 
of the two source classes fluctuate, however, this quan- 
tity can go from positive to negative, forcing us to switch 
the choice of signs when this happens. Notably, a sign 
change would have to occur when /iC^^ = l2Ci.2, pre- 
cisely when C^.tot is minimized. Indeed, since I2 / Itot is 
changing through the minimum, the sign must change 
at every minimum with the only exception being when 
I2/ hot is simultaneously minimized or maximized, which 
does not happen in general. Consequently, no single 
choice of signs describes the entire spectrum, but only 
the region between two consecutive minima. 
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F. Decompositions from Multiple I s 

In the case where two distinct anisotropy energy spec- 
tra, each measured at a different I, can be obtained, a 
full decomposition is possible provided Ci varies with I 
for at least one source class and can be determined 
for each spectrum, e.g. via a plateau. We still assume 
Ci is independent of energy. Differentiating Eq. [3] with 
respect to energy and rearranging, we obtain. 



dE [hot J 2{^[Ce.j{l) + CiM ~ Ci,2m 



(25) 



Since the left-hand side is independent of i, we must have 
for two different £ values £i and £2 at any energy E 



dCi,^tot/dE 



(26) 



rfCf2,tot/ dE 
2[-£7(C^2,i + C'^2,2) - C^2,2] 

Using Eq. [3]to eliminate €^^^2 and Ct^,2 and solving 
for Ii yields 



dC, 



I2 .tot 



dE 



a 



dCt^_ 



^2, tot 



dCi2 



dE 



Ce. 



(27) 



dE 



This method is particularly ideal in that the input 
spectra are derived exactly and without making any as- 
sumptions about their relative intensities or anisotropics. 
However, this method may be difficult to implement for 
the IGRB because we expect most gamma-ray source 
populations to produce fluctuation angular power spec- 
tra that are dominated by the Poisson angular power, 
which takes the same value at all £ (e.g., 25, 31, 32, 58]), 
with the notable exce ptio n of dark matter annihilation 
or decay (e.g., [H [H M HI]). 



V. THREE COMPONENT DECOMPOSITION 

The separability of the equations for the total inten- 
sity and anisotropy at a given energy into contributions 
from each component source class enables a background 
of multiple components to be decomposed provided all 
source classes are uncorrelated and the contributions of 
all but two source classes are known. Furthermore, un- 
der certain conditions a decomposition can be performed 
even when information about additional source classes is 
more limited. 

In this section, we discuss specifically the case of a 
three-component background. The intensity of a three 



component background as a function of energy is given 
by 

hot ^h+h + h (28) 
and the fluctuation anisotropy as a function of energy by 

Ct,tot{E)^(-^\ Ci.i+(^) Ci.2+(^) Ci,z 

V^tot/ V^tot/ \hot J 

(29) 

again assuming uncorrelated components. 

In the following, we examine the applicability of the 
techniques we discussed in two plausible scenarios for the 
IGRB composition. 



A. One component with known intensity and 
anisotropy 

It is possible that we can obtain, via some other analy- 
sis, expressions for both the intensity and anisotropy of a 
third component. Is and (7^,3 as functions of energy. Such 
a situation could occur for a component arising from a 
population of (almost) uncorrelated bright point sources, 
for which enough individual members have been resolved 
so as to obtain a thorough understanding of the spectral 
behavior, and to constrain well the distribution of source 
fluxes, dN/dF, down to a point from which a reasonable 
extrapolation to even lower fluxes is possible. One can 
envision, for example, this to be the situation for gamma- 
ray blazars after the completion of the Fermi mission. 
From dN/dF both the anisotropy level and the overall 
intensity normalization can be calculated, and from the 
understanding of individual source spectra the energy de- 
pendence of the intensity can be evaluated. 

We can thus rewrite our equations as 



(30) 



and 



e.tot 



h 



a 



e.2 ■ 



(31) 

where the left hand side of both equations are deter- 
minable directly from observables and knowledge of the 
third component's properties. Hence, we have reduced 
the problem back to the two-component case where all 
of our decomposition methods apply. 



B. One component with zero anisotropy, and a 
second component with known spectral shape 

A second likely three-component scenario that is work- 
able in this formalism is one in which no component is 
completely known, but where one component (compo- 
nent 1) has a well- understood and zero (or negligible) 
anisotropy, and a second component (component 2) has a 
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well-understood intensity spectral shape, even if its over- 
all intensity normalization is unknown. 

In the context of the IGRB, the zero-anisotropy com- 
ponent could be contamination from unrejected cosmic- 
ray electrons entering the detector [20} or a combina- 
tion of such cosmic-ray contamination and a very-low 
anisotropy cosmic component (e.g., star- forming galax- 
ies, or cascade emission in the case of significant inter- 
galactic magnetic field [60]). Blazars, on the other hand, 
could be the component with a well-understood spectral 
shape (since, by the end of the Fermi mission, thousands 
of blazars will have been resolved and have their individ- 
ual spectra measured) , even if the overall normalization 
of their intensity contribution to the IGRB is still uncer- 
tain. 

We will see that in this case the energy dependence of 
the intensity of a third component can be determined up 
to a normalization constant. Such information could have 
extremely high impact if the third component (compo- 
nent 3) is, for example, a contribution from dark matter 
annihilation or decay, as we discuss in WII 

Because component 1 has zero anisotropy, Eq. (|31|) be- 



h{E) 
hot{E) 



IsjE) 
hot{E) 



(32) 



where we have written explicitly all energy dependencies. 

We assume, as before, that there is an energy where 
we know that component 3 does not contribute signifi- 
cantly (as could be the case at low energies for certain 
dark matter annihilation or decay components). At i?o, 
Eq. ([5^ then becomes 



(33) 



Now since component 2 is assumed to have a known spec- 
tral shape, we can write 



hiE) = l2iEo)giE) 



(34) 



where g{E) is a known function of energy, normalized so 
that g{Eo) = 1. Solving Eq. §^ for I^iE) and using 
Eqs. (|33| and ()34|) to eliminate Ci^2 and hiE), we then 
obtain 



ll{E)C,,:, 



{E)Ce.tot{E) - g\E)Cttot{Eo)lLiEo) ■ 

(35) 



Since all quantities on the right hand side of Eq. ([35]) 
are known, it follows that I^{E) can be derived up to a 

normalization constant l/\/ Ci,^'- 



h{E) 



I lldE)Ci,tot{E) 



g^E)Ce,tot{Eo)lLiEo) 



C, 



(36) 



VI. DISCUSSION 

Unraveling the contributions of multiple source classes 
to a diffuse background can be accomplished in many 
cases by combining intensity and anisotropy informa- 
tion. We have presented model-independent decompo- 
sition techniques which can recover the spectral shapes 
of the constituents of a two-component diffuse back- 
ground (low-anisotropy plateau, known zero-anisotropy 
component) and techniques which can recover both the 
shape and normalization of the component spectra (dou- 
ble plateau, high-anisotropy plateau, minimum, multi- 
ple t). Additionally, we have discussed cases in which 
these techniques can be applied to backgrounds of three 
or more components. 

The techniques presented here are applicable at any 
wavelength and for any diffuse background composed 
of uncorrelated source classes (for example, in the case 
of the microwave background, the cosmic early-universe 
component, the dust emission from star- forming galax- 
ies, and the contribution from blazars are uncorrelated 
or very weakly correlated components). Because each 
method has a clear signature in the observed anisotropy 
energy spectrum and/or straightforward mathematical 
tests of validity, there is no ambiguity in selecting which 
technique to apply. With the exception of the known 
zero-anisotropy component technique, no a priori as- 
sumptions need to be made about the nature of the com- 
ponent spectra. Hence, the methods presented here are 
unique and mo del- independent. 

Each technique detailed in this paper, however, re- 
quires the source classes to have energy-independent fluc- 
tuation angular power, which is not a perfect assump- 
tion for many astrophysical sources, extragalactic ones 
in particular. Such a dependence would introduce de- 
generacy into the decompositions as any variations in 
the anisotropy energy spectrum could be explained ei- 
ther by changes in the fractional makeup of the intensity, 
as we have assumed here, or by changes in Ci with E 
for each source class. However, it is likely that these as- 
sumptions hold over at least some energy ranges. While 
in certain cases these methods can be applied to scenar- 
ios in which more than two components contribute to 
the diffuse emission, not all scenarios with three or more 
relevant contributors can be decomposed using the tech- 
niques presented here. Thus, a significant contribution 
from a third component in scenarios other than the two 
cases discussed in f|V]may render the application of these 
methods difficult at some energies. 

The ability to decompose the intensity energy spec- 
trum of the IGRB and recover the constituent spectra, 
even at an accuracy of up to a multiplicative constant, is 
of paramount importance in understanding the physical 
properties of the underlying source classes. Specifically 
in the case of blazars, the slope of their gamma-ray back- 
ground contribution reveals the^spectral properties of the 
unresolved blazar population [36] or the relative contribu- 
tion of different type of blazars to gamma-ray background 
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and consequentl y to the faint end of the blazar luminos- 
ity function 'i^, [S^l ■ If blazars are indeed a subdomi- 
nant component of the gamma-ray background intensity, 
as suggested by recent constraints from the measured 
IGRB anisotropy iS^, such a decomposition would 
in principle allow us to deduce the slope of their collec- 
tive emission at much higher accuracy than by modeling 
and subtracting the dominant components, the details 
of which may be largely unknown. In addition, a decom- 
posed intensity spectrum for blazars can place constraints 
on the intensity of the EBL. In the case of dark matter, 
a decomposed intensity spectrum is an uncontaminated 
measurement of the photon spectrum from dark matter 
annihilation or decay, which in turn can provide informa- 
tion about the dark matter particle mass and dominant 
annihilation or decay channels. 

At the same time, several of our techniques can con- 
strain the source population anisotropy as well. Such 
constraints are also extremely important in understand- 
ing the statistical properties of a source class, as they 
provide information about the faint end of the luminos- 
ity function that is independent from that encoded by 
the collective intensity from unresolved members of the 
class. 

Taken in complement with other analysis methods, 
these techniques can provide a unique, valuable window 
through which to probe the physics of the IGRB or any 



other diffuse astrophysical background. 
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FIG. 4: Example known zero-anisotropy component decom- 
position. Top: Total IGRB intensity (top) and anisotropy 
(bottom) energy spectra as measured by Fermi-hKT using 
~ 21 months of data. If the data point is within 3cr of zero, 
we place a 3a upper limit bar in addition to the la error 
bars. Since the normalizations of the intensity spectra are 
not recovered with this technique, the spectra are shown with 
arbitrary normalization. Bottom: The decomposed intensity 
energy spectra of the anisotropic component (blue circles) and 
the zero-anisotropy component (red triangles) are recovered 
using the known zero-anisotropy component technique. Each 
component's input intensity spectrum is overlaid in black. 



